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Abstract. The abundance and distribution of collapsed objects such as galaxy 
clusters will become an important tool to investigate the nature of dark energy and 
dark matter. Number counts of very massive objects are sensitive not only to the 
equation of state of dark energy, which parametrizes the smooth component of its 
pressure, but also to the sound speed of dark energy as well, which determines the 
amount of pressure in inhomogeneous and collapsed structures. Since the evolution of 
these structures must be followed well into the nonlinear regime, and a fully relativistic 
framework for this regime does not exist yet, we compare two approximate schemes: the 
widely used spherical collapse model, and the pseudo-Newtonian approach. We show 
that both approximation schemes convey identical equations for the density contrast, 
when the pressure perturbation of dark energy is parametrized in terms of an effective 
sound speed. We also make a comparison of these approximate approaches to general 
relativity in the linearized regime, which lends some support to the approximations. 
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1. Introduction 

We now have overwhelming evidence that the Universe is accelerating, possibly under 
the influence of some type of negative-pressure substance - dark energy (DE) [HE], [3]. 
However, even though DE may be directly responsible for this enhanced expansion, it 
is widely believed that the direct impact of perturbations in DE density and pressure 
on structure formation is very weak. This is strictly correct only for a cosmological 
constant model of DE, which does not have perturbations. 

For most scalar field models of DE, this component remains very homogeneous even 
on galaxy and cluster scales. Heuristically, this can be understood as follows. In these 
models the scalar field can not have relaxed to its minimum energy state and one must 
require that the time scale for the variation of the field is longer than the Hubble time, 
implying a very flat potential. Therefore, the scalar field must be extraordinarily light, 
m < H 0l where H is the Hubble parameter today. The mass of the scalar field sets the 
scale for its spatial variation and hence one usually expects small perturbations in the 
scalar field for scales A < 1/m (the Compton wavelength) , which are of the order of the 
Hubble radius. However, this argument may not apply to more general models of dark 
energy. 

If our only concern is the evolution of the background, then the role of dark energy 
in the evolution of dark matter perturbations is completely determined by its equation of 
state w = p e / p e , where p e is the homogeneous pressure and p e is the homogeneous energy 
density of dark energy pO, [5J E]. At this level, dark energy affects structure formation 
indirectly because, as it starts to dominate the background, very large structures are 
ripped apart by the ensuing accelerated expansion [2 [8]. 

However, dark energy can influence structure formation in an additional manner. 
If it is a dynamical field or fluid, then dark energy must possess inhomogeneities, and 
these perturbations will interact gravitationally both with themselves and with clumps 
of dark matter [9j. This means that, unless dark energy is just a cosmological constant, 
it will both feel and create local gravitational potentials. 

Although the effect of these inhomogeneities in the dark energy component becomes 
small as w — > — 1, in many models with w ^ — 1 it can be non-negligible when evolved 
in the nonlinear regime P3H [HJ H21 UBl U31 UBl US1 HH UB1 USl 1201 EH 1221 E3J - Since 
the effects of dark energy perturbations on the cosmic microwave background are quite 
small (see, e.g. [21]), structure formation is the only remaining probe of the nature of 
dark energy on small and intermediate scales. 

Nevertheless, a fully relativistic method to treat nonlinear perturbations is not 
available. When there is a pressure ingredient the nonlinear relativistic equations take 
a very complicated form. The Lemaitre-Tolman-Bondi (LTB) model [25], [26], [27] is the 
closest one can get to a working formalism, but it only works if matter is pressureless. 
The problem is not with the gravity side of the equations, but with the nonlinear 
evolution of matter and the relativistic treatment of pressure. 

In this respect, the only well-studied models with inhomogeneous dark energy are 
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those involving canonical scalar fields [H3, [HI [TZI [THl [D2, ED], for which the equations 
of motion and the pressure follow directly from a given Lagrangian. For these models 
the free parameters are the scalar potential and some set of initial conditions. In this 
approach, the equation of state, the density perturbations and the pressure perturbations 
are derived quantities. Hence, a more kinematical and model-independent approach to 
structure formation, closer in spirit to the homogeneous description of dark energy in 
terms of a parametrized equation of state w(z), is sorely lacking. 

There are two very different approximations to full-blown general relativity that 
have been frequently used. They are the spherical collapse (SC) model [28, 29, 30j and 
the pseudo-Newtonian (PN) approach [3H [32], [33j [34]. We have recently used these 
approximations in the nonlinear regime in order to show that the "effective equation of 
state" of dark energy inside a collapsed region could be very different from its background 
value [23] . 

In this work we show that, even though the underlying assumptions for either 
approach are rather different, they yield exactly the same nonperturbative equations 
as long as the pressure perturbations are treated in the same way. They also have 
an important advantage: they allow for a completely parametrized approach to dark 
energy. Furthermore, we compare the growth of perturbations in the linear regime with 
a linearized relativistic analysis and show that they are similar, lending support to the 
approximations. 

This paper is organized as follows. In Section 2 we review both the PN and SC 
approaches and show that they are equivalent. In Section 3 we study the linear evolution 
of perturbations in DE in this approximation. The linear evolution of perturbations in 
a universe with a 2 component fluid is studied in general relativity in Section 4. We 
present a comparison between the relativistic analysis and the approximate analysis in 
the linear regime in Section 5. Section 6 concludes. 

2. Spherical collapse and pseudo-Newtonian cosmology 

In linear perturbation theory there are essentially three degrees of freedom for scalar 
perturbations: the energy density perturbation 5p, the pressure perturbation 5p and the 
scalar anisotropic stress n [35, 36j. An alternative set is given by the density contrast 
Sp/p, the velocity potential 6 = V ■ v and the anisotropic stress [37]. Since large-scale 
anisotropic stresses decay rapidly, they can only become relevant again inside structures 
which have collapsed. This means that anisotropic stress should not influence the mass of 
these structures, and therefore it is unlikely that dark energy models can be differentiated 
on the basis of anisotropic stress. For this reason we do not consider it any further in 
this work (see, however, [38J). 

We will parametrize the pressure perturbation using the so-called effective sound 
velocity [39], defined as Cg ff = Sp e /Sp e . We will assume that Cg ff is a function of time 
only, even though this simplification lacks any formal basis in cosmological perturbation 
theory. This should be clear from the fact that 5p e is an independent degree of freedom 
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whose time and spatial dependences can be, and often are, completely different from 
5p e . Only in a particular gauge (the so-called "rest frame" of the fluid, where Tq = 0) 
the effective sound speed coincides with the universal sound speed of linear relativistic 
perturbations, c 2 x [39~ll40|. It may be difficult to realize this parametrization in a natural 
model, but the situation is not much different from what happens when we parametrize 
the equation of state. 

The main reason that we use the effective sound speed, though, is that it allows us 
to study nonlinear structure formation within the spherical collapse model [28J. In this 
extremely simple model, a spherically symmetric region of homogeneous overdensity 
evolves inside the homogeneous expanding Universe (this is the so-called "top-hat" 
density profile). General relativistic arguments show that one can regard the overdense 
region as a mini-universe of positive curvature, and then we use the Friedmann and the 
Raychaudhury equations to evolve the density and radius of the spherical region [29, 30J. 

It is therefore extremely interesting that this simplified relativistic approach 
coincides with a pseudo-Newtonian approach to cosmology. In fact, we will show below 
that, as long as the pressure perturbations are described in terms of an effective pressure, 
the two approximations are completely equivalent. This means that the main physical 
characteristics of gravitational collapse of structures such as galaxy clusters is probably 
well described within this framework. 

The argument is as follows. First, the SC approach should be a good approximation 
for large scales (where relativistic effects should matter most), but not necessarily for 
small scales, where the "mini-universe" argument is less persuasive. On the other hand, 
the PN approach is well-motivated by the physics of gravity in small scales, but is not 
assured to work for large scales. That the two approaches coincide shows that, at least in 
some limited sense, the equations of the SC/PN approach should give a good description 
of the gravitational interactions on scales smaller than the Hubble radius. 

2.1. Pseudo-newtonian cosmology 

In PN cosmology, particles in a comoving grid attract each other gravitationally with 
a Newtonian potential. The positions of the particles in the grid are the perturbed 
variables. Although obviously limited, this approach can be used for any configuration, 
not only the spherically symmetric ones. But in order to bring the PN approach closer 
to the SC model, we will adopt the same basic assumptions of the SC model for the PN 
cosmological perturbations. 

We consider an admixture of two fluids, cold dark matter and dark energy. The key 
assumptions of the SC model (see the next subsection) are that the density of each fluid 
is homogeneous at all times in the spherical region (this is the top-hat density profile), 
and that the velocity profile preserves this homogeneity. 

The comoving coordinates are xq = fo/a, where r*o is the homogeneous 
(unperturbed) physical distance - here, the radius of a spherically symmetric region. 
Under the assumption of the SC model, the perturbed physical distance (physical radius) 
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can be written as: 

f=[a{t) + f{t,x )]x , (1) 

where a is the usual scale factor and / is the function that accounts for the deviations 
from the background evolution. The physical velocity is then given by: 

where ' = d/dt and H = a/a is the Hubble parameter. From the last equality we can 
define an effective rate of expansion for the spherical region: 

h = H + L (3) 
a 

Since the perturbed velocity is related to the peculiar velocity v by 

u = a x + v , (4) 

we obtain from Eq. ((2)) that: 

v = fx . (5) 

In particular, the divergence of this velocity field is given by: 

e = V-v = 3f + x -Vf . (6) 

But for a top-hat profile the last term vanishes, and we obtain a simple relation between 
the local expansion rate h and the background expansion rate H: 

h = H + t = H+-H- . (7) 

a da 



The PN cosmological model is described by the equations [31] : 

?^+V-(u j p j )+p j V-u j = 0, (8) 



a + ^- v J^=- v *-^' (9) 



V 2 $ = 47rG^(p fc + 3p fc ) , (10) 

k 

where pj, pj and Uj denote, respectively, the density, pressure, velocity of a given cosmic 
fluid and $ is the Newtonian gravitational potential due to all the components; the 
equations are written in physical coordinates. The corresponding perturbations above 
the background are denoted by 8pj, 5pj, Vj and 0. These equations are, respectively, 
generalizations for fluids with pressure of the continuity equation, of the Euler equation 
(both valid for each fluid species j), and of the Poisson equation. Notice the absence of 
an equation that dictates the evolution of pressure: in this hydrodynamical approach, 
pressure is a thermodynamical function of the energy, temperature, etc. 

For cold dark matter and baryons the pressure is zero, but for dark energy there 
is a homogeneous as well as an inhomogeneous pressure. The homogeneous pressure 
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is usually described in terms of a parametrized equation of state w e (t), such that 
p e (t) = w e (t)p e (t). As for the pressure perturbations, we have chosen to specify another 
free function, the effective sound speed Cg ff , so 5p e = c^. s 5p e . Within the SC description, 
this means that we consider an effective equation of state w c inside the spherical region 
which is not necessarily equal to the background equation of state. 

With the assumptions of the SC model, the equations of PN cosmology assume a 
simple form. Using the density contrast 5j = d~Pj/pj we obtain, after some algebra: 

5j + 3H (c c 2 ff j - Wj ) 5 3 + e -i[l + Wj + (l + Ce 2 ff ,) Sj] = , (11) 

CI 

Q2 

6 3 + H9j + J- = -AirGa ^ p k 5 k (l + 3c c 2 ff fc ) . (12) 

k 

Eq. ( TlTi ) follows from the continuity equation, and Eq. ( IT2l is the divergence of the 
Euler equation. The last equality in Eq. (fT2l is found by using the Poisson equation. 
Note that, in general, we have separate Euler equations for each fluid [21], but for a 
top-hat profile (V<5j = 0) they turn out to be identical, so there is only one 6. The 
reason for that is obvious: in order to preserve the top-hat profile, all fluids must flow in 
the same way. Hence, in this approximation we have something similar to an effective 
single fluid description |41j. 



2.2. The spherical collapse model 

Let us now briefly review the spherical collapse model. This formalism describes a 
spherically symmetric region of uniform energy density p c = Po + $P immersed in a 
homogeneous universe of energy density po- This spherical region will detach from the 
expansion of the Universe and eventually collapse. 

Consider the continuity equation for each fluid denoted by an index j in the spherical 
region: 

p Cj + 3h (1 + w Cj ) p Cj = , (13) 

where h — r/r is the local expansion rate of that region and w Cj denotes the equation 
of state in the perturbed region. We can regard this spherical region as a Friedmann 
Universe with spatial curvature |28| . The dynamics of the coordinate r is then given by 
the second Friedmann equation applied to this collapsing region: 

- = "— E(^ + 3 ^) • (14) 



Equations (1131 ) and ( fHj) , which were obtained using general relativistic arguments, are 
the basic equations of the SC model. Note that there is only one dynamical equation 
for the collapsing region, which is in agreement with the single Euler equation that we 
found for the velocity field in the PN description, Eq. (fT2l . 

The pressure and the energy density outside the spherical region are related by the 
background equation of state, po j = w 0i p . Inside the spherical region these quantities 
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can be different from their background values, so we have now p Cj = w c .p Cj for the 
collapsing region. In order to compare the SC formalism with the PN equations derived 
in the last section, we will employ here the same effective sound speed we used before 
in order to describe the pressure perturbations. Hence, we need to express the equation 
of state w Cj in terms of Cg ff j. Using the density contrast Sj = Spj/po^, we have that: 

Pc j = (l + 5 j )p Qj , (15) 

from which it follows that: 

p C] _ Po, + S Pj _ - W .)-*L- (i 6 ) 

This equation relates the equation of state in the perturbed region to the background 
equation of state, the effective sound speed and the size of perturbations. It is possible 
that the nature of dark energy can be significantly changed in the perturbed region, a 
phenomenon we dubbed "dark energy mutation" [23J. This effect is general, occurring 
even at the level of linear evolution, and its magnitude depends on the dynamical 
evolution of DM and DE fluctuations - see also Refs. [TBI 120]. 
Using now Eqs. (fT5|) and ([TBI) we can recast Eq. (fl3|) as: 

6 . + ( 3/i _ 3H) (i + w .) (i + 5 .) + 3h ( c 2 ff . _ Wj ) Sj = o. 
We can eliminate h using Eq. fTJ), with the result: 

^ + 3H {c 2 cS] - Wj ) 5, + [1 + Wj + (1 + c e 2 ff Sj] - = . (17) 

(X 

Now consider the dynamical equation (fT4l) . From Eq. j7|) we can write: 

h = ¥ --h* = H+±- -H(- , (18) 
r 3a oa 

and substituting this expression into Eq. (|T4l) we obtain, with the help of Eqs. (fT5ll -(fT6ll. 

that: 

9 + H9+^ = -AnGaJ^PoJk (l + 3c c 2 ff k ) . (19) 

k 



Equations (1171) and ( [1911 are identical to Eqs. (jTTi)-(jT2l). This means that both 
approaches are identical. The relations and (|T6l) enable us to translate the PN 
variables into the SC variables, and now it becomes clear that the two different 
descriptions give the same dynamics for a top-hat perturbation where pressure gradients 
are absent. 



3. Linear evolution in the SC/PN approach 

Even though we showed that the PN and SC approaches are equivalent, that still does 
not mean that they are correct. Unfortunately, presently there is no fully nonlinear 
general treatment of the evolution of perturbations in General Relativity (GR). For this 
reason, we will compare our linearized results with those obtained from linearized GR. 
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We will compute the linear evolution of an overdense region well inside the matter- 
dominated era, and will compare the growing mode obtained in the PN/SC formalism 
with the relativistic growing mode. 

The first-order equations can be linearized and recast as a single, second-order 
differential equation for the density contrast of each fluid species. We will assume that 
there is always a dominant (d) and a subdominant (s) fluid. Using the scale factor a for 
the time evolution (' = d/da), we obtain for the dominant species: 



3A d + - (1 - w d ) 



(20) 



3w d )-(l + w d )(l + 3c 2 cSd )}=0. 



where 



(21) 
(22) 



A d = (Ccff d ~ W d ) . 

For cold dark matter (c e g — w — 0) this equation reduces to 
3^ 35, 
6d + 2 a + 2«2 " U ' 

with the well-known growing solution 5(a) oc a. Hence, when cold dark matter is 
dominant, which should be the case in the linear regime, the linear evolution of its 
density perturbations is the standard one. 

For the more general case, of a dominant fluid with constant equation of state and 
constant speed of sound, the solution is given by: 



5 d (a) = Cl a 1+3wd + c 2 a- 3 ( 1+2c ^ 



-w d )/2 



(23) 



Turning now to the the sub-dominant fluid, its perturbations obey the equation: 



a 



3A s + -(l-w d ) 



(24) 



35, 
"2a 2 



[A,(l-3^)] = ^(l+ Ws )(l + 3c c 2 ff J , 



where 



eff s 



W 



) • 



Assuming again that w s and c^ Ss are constants one has the solution: 



$s (a) = c 3 a 



-3A S 



C4a (Sw s -l)/2 + C5a l+3 Wd 



where 



c 5 = 77 c i (! + w s) (1 + 3c 2 ffd ) 



(25) 
(26) 
(27) 



x 



: I 3«v) ( A 3 + i(l + w d U + X - (1 - 3w d ) A s 



~i -i 



arises from a particular solution of the inhomogeneous equation. 
In particular, if matter is the dominant fluid it follows that: 

ci (1 + w s ) 



c 5 



3A S + 1 



(28) 
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and the dark energy density contrast grows in the same way as the dark matter density 
contrast. In addition, for the case in which c 2 s = w one has an adiabatic condition 
satisfied, namely 5 e — (1 + w e )5 m . In general, however, the perturbations have a non- 
adiabatic component and the dark energy density contrast evolves as: 

5 S (a) = Tj^qpyfc (°) + c 3«" 3As > ( 29 ) 

where the last term in the right-hand-side is a decreasing mode in most cases. 



4. Linear Evolution in GR 



In a previous paper |21| we showed that, for a single perfect fluid with no pressure 
gradients, the growing modes in the linearized SC/PN approach coincide with those 
found with General Relativity (GR). Now we want to compare the PN and the GR 
solutions for the dark energy perturbations in the linear regime, including pressure 
gradients. We will consider these perturbations during the matter-dominated period, 
i.e., while DE is subdominant. This is motivated by the fact that most observed 
structures were formed well into the matter-dominated period, 

We consider scalar perturbations to the metric in the newtonian gauge without 
anisotropic stress: 

ds 2 = (1 + 2(p)dt 2 - a 2 (1 - 2(p) dx 2 . (30) 
The (00) and (ii) componentes of Einstein equations in Fourier space are: 

^<f> + 3HU + H<f)) = -4irGj2 5 Pj > (31) 

(X \ / 

3 

+ + (2- + H 2 J <f> = 4:irGj2 6 Pj > ( 32 ) 

3 

and the conservation equations Tq.^ = and = yields: 

63 + 3H (c4 i - Wj ) Sj + (1 + Wj ) - 3</>J = , (33) 

+H(1- 3cM 6j - — ^ -0 = 0, (34) 

3 v S3 ' 3 (1+Wj)pja a v ' 

where c 2 SJ - = pj/pj is the adiabatic speed of sound. . 

In summary, the evolution of perturbations in a system consisting of dark energy 
and dark matter in linearized GR is described by the following set of 5 coupled 
differential equations: 

+ 4#0 + (2°- + H 2 ^<P= ^H 2 n e c 2 e{i S e , (35) 

L + — - 30 = , (36) 
a 
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5 e + (l + w e ) - 3<fij + 3H (c c 2 ff - w e ) 6 e = , (37) 
k 2 

6 m + H6 m = 0, (38) 

a 

e e + H(l-3c 2 se )9 e - .ff* 5 ? --0 = 0. (39) 
' [l + w e )a a 

5. Comparison between GR and PN 

In PN cosmology the linear evolution of DM and DE is determined by the system of 
equations that arise from Eqs. (JSj) and ([9|) for each fluid, namely: 

L + — = , (40) 
a 

S e +(l + We) -+3H (C 2 CS -W e )S e = 0, (41) 
(X 

k 2 

6 m + H6 m -—<j> = 0, (42) 
a 

k 2 r 2 k 2 

9e + H6 e - *y d ° --0 = 0. (43) 
(1 + w e ) a a 

Notice the absence of a dynamical equation for 0. To eliminate the k 2 (p terms 
we can use the constraint implied by the Poisson equation in PN cosmology, Eq. ffTOl . 
Then the time variation of the potential is determined by the evolution of the density 
constrasts. Also notice that these equations lack some terms when compared with their 
relativistic counterparts, as already pointed out in Ref. [42]. However, as we will show, 
during the matter-dominated era and on small scales, this discrepancy changes only the 
velocity that DE perturbations decay but do not modify its late-time behaviour. 

In the matter-dominated regime, = const, is a solution of Eq. (l32l) . which also 
arises from the system of Eqs. ([40l - (l43l) . In this case, it is interesting to notice that 
Eq. (j37l) becomes identical to Eq. (l411 . However, Eq. (|39l coincides with Eq. (T43l only in 
the case c se = 0. 

As we see, the equations for the growth of perturbations are different in GR and 
PN already in the linear regime. Now we perform a quantitative study of this difference. 
We will work out the case of a matter-dominated universe with a small DE component, 
as expected in the linear regime, in which case is a constant. Furthermore, to avoid 
further complications, we assume constant values for w and Cg ff . 

Under these conditions we can write a second order differential equation for the 
linear growth of the dark energy density perturbation as a function of the scale factor, 
5e(a): 

6, 

a 2 H 2 a" a" a" a 



5' 
a 



p 2 



v 1 + w )-r7^- ( 44 ) 
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This equation arises both in GR and PN: in the latter case, we keep the pressure gradient 
in the Euler equation ([9]), which was dropped in the case of a top-hat perturbation. Only 
the parameters a and f3 are different in the two cases: 

a GR =^ + 3A-3w- /3 GR = 3A Q - 3w J (45) 
a PN = g + 3A ; (3 PN = ^A , (46) 



where A was defined in Eq. (i25l ). 

We will make a comparison focusing on small scales, where the PN approximation 
is supposed to be more accurate. In this case, we can neglect the (3 term in the square 
brackets of Eq. (l44l) . and we immediately write a constant particular solution: 

Se = (47) 

C eff 

In order to solve the homogeneous equation we perform the following change of 
variables: 

5 e (a) = x 1 - a y(x) , (48) 
where x is defined in terms of the conformal time r\ as x = kc e ^r]. Then Eq. (T44l) becomes: 



d 2 y 1 dy r 



dx 2 x dx ^ 



x' 2 



y = , (49) 



where, according to the different coefficients in Eqs. (f4"5|) - (l4l)]) . \i assumes different 
values: 

^ = ±-(1-64) , (50) 

fi PN = ±- (l + 6w - 6c 2 cS ) . (51) 

The solutions are Bessel functions of first kind J± M (x) . The dark energy density 
contrast behaves as: 

6 e (x) = x 1 ^^, (x) - 0. (52) 

C eS 

These solutions both have an oscillatory behaviour with a decreasing amplitude 
proportional to x 1_a_1 / 2 and they eventually reach the constant value S e = 
-(l + W )0/c e 2 ff . 

In order to check this analytical behaviour we numerically solve the complete system 
of coupled differential equations (I35l)-(i39l. We used as illustration c 2 s = —w e = 0.8, 
Qjfe = 1 — £lm = 0.75, and we evolved the equations from an initial redshift of = 100. 
We examined the mode k = 100H = 0.0236/iMpc -1 , corresponding to a physical scale 
of A = 266/i~ 1 Mpc, well inside Hubble radius at Z{ and large enough to be in the linear 
part of the matter power spectrum. As initial conditions we chose 0, = — 10~ 4 , fa = 
and: 

kHl + z,) 2 - 



1+ 3H{ Zl ) 2 



8 e {zi) = (1 + w)5 m (zi) ; (53) 
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, (54) 



2(1 + Zi )k\ 
Om{Zi) = (pi ; e e (Zi) 



which are consistent with Einstein's equation and adiabaticity. The result is presented 
in Fig. (1) and compared to the decay factor and the final value given by Eq. fl47l) . 



o 




0.25 



Figure 1. Linear evolution of dark energy perturbations in GR at small scales 
(k = 0.0236/iMpc -1 ) for c 2 eS = -w e = 0.8. The solid line is the solution of the 
complete set of 5 coupled differential equations. The dotted line is the decay factor 
according to Eq. (]52]) . The dashed line is the particular solution Eq. ([47]) . 



We also perform the same exercise for the PN approximation. The numerical 
solution of the system of Eqs. (I40~l) -(l43l) with the same parameters and initial conditions 
is presented in Fig. (2) and compared to the decay factor and the final value given by 
Eq. fliTI) . Again we see that the qualitative analytical behaviour is reproduced by the 
numerical solution. 

Therefore, even though the equations from GR and PN are not the same already at 
the linear level, the results are not qualitatively different. In particular, both approaches 
predict the same asymptotic behaviour for the DE perturbation. Because the decay rate 
of the transient is slightly different in each case, the time when the asymptotic regime 
is reached differs - in the PN approach this happens at a later time. 

At this point we should call attention to the origin of an apparent discrepancy 
between the results obtained in this Section, namely, a constant behaviour of the DE 
perturbations, Eq. ([47]), and the result obtained in Section 3, where we showed that 
in the SC/PN approach the DE perturbations grow as DM perturbations, Eq. ([29]) . 
The reason is that in Section 3 we assumed a top-hat profile of the perturbation, which 
amounts to setting kc e ^ = in the square brackets of Eq. ([44]) . In this case, the 
particular solution is 

l + w k 2 <f) 3(1 + w) 



5 e (a) 



(3 H 2 a 2 



2(3 



(55) 
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Figure 2. Linear evolution of dark energy perturbations in the PN approximation 
at small scales k = 0.0236/iMpc _1 for c 2 s = —w e = 0.8. The solid line is the 
solution of the complete set of 4 coupled differential equations. The dotted line is 
the decay factor according to Eq. (}52]) . The dashed line is the particular solution 
Eq. (j47j). 



where in the last equality we used Poisson's equation for the case of a dominant dark 
matter component, k 2 <p = — (3/2)H 2 a 2 5 m . Hence, we see that indeed in this case, or in 
fact for perturbations with a small mode number k, the perturbations in DE grow at 
the same pace as the DM perturbations in the linear regime. 

On super-Hubble scales, PN cosmology is not expected to be valid, due to its 
inherently instanteneous interactions: indeed, in that framework perturbations with 
scales larger than the Hubble radius would behave in the same way as those well inside it. 
However, since we are only interested in the evolution of perturbation which are initially 
in the linear regime and well inside the Hubble radius, this mismatch is irrelevant. 
Therefore, we do not compare the PN perturbations with the GR perturbations in large 
scales. 

As a final remark, we recall that the analytic solution, Eq.(|52l). is valid only for 
linear perturbations during the matter-dominated period. In this regime the matter 
density constrast grows as 5 m oc a and is constant in time. When the structures enter 
the nonlinear regime, matter fluctuations must grow faster, i.e, 8 m oc a n , with n > 1, 
then the gravitational potential should grow in time. Hence DE behaviour in nonlinear 
structures is expected to be different from the linear analysis results. However, the 
asymptotic constant solution in Eq. (l52l) is valid during the initial nonlinear process of 
matter collapse and DE fluctuations can grow with 0. 

6. Conclusions 

The study of perturbations in dark energy has received a great deal of attention 
recently. DE perturbations have the potential to alter the process of large scale structure 
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formation in the universe. The existence of DE perturbations can in principle be tested 
in future surveys such as the Dark Energy Survey (DES) [43] and EUCLID |44j. These 
future observations may help to distinguish among different models of DE. 

Structure formation occurs during the nonlinear stages of the evolution of 
perturbations. Unfortunately, there is no rigorous analytical description of this nonlinear 
stage in full GR. Aproximation methods must be used. Possibly the most trusted method 
is N-body simulations, but due to its very intensive computing requirements, it is not 
practical when one wants to study different models. Furthermore, N-body simulations 
employ newtonian physics and do not allow for the possibility of DE fluctuations. 

In this paper we study two different approximation schemes, namely the Spherical 
Collapse and Pseudo-Newtonian approaches. The advantage of these schemes is that 
DE can be fully characterized by 2 functions: the equation-of-state parameter w(z) and 
the effective speed of sound c e g(z). We show that, under a minimal set of assumptions, 
it is possible to translate one approach into the other, rendering them completely 
equivalent. In order to compare these approximations with GR, we study perturbations 
in the linearized regime with all approaches. When the assumptions about the pressure 
perturbations are the same both in GR and PN/SC we find that the fluctuations present 
the same qualitative behaviour, lending support to the approximations. However, in 
order to establish more firmly the validity of the approximations in the nonlinear regime, 
a comparison should be made with some nonperturbative model in GR, such as an 
extended LTB class of models, including fluids with pressure. Work along this direction 
is in progress. 
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